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Abstract 

The parity nonconservation (PNC) effect on the laser-induced 2 3 S*i- 2 1 S'o transition in heavy heliumlike 
ions is considered. A simple analytical formula for the PNC correction to the cross section is derived for 
the case, when the opposite-parity 2 1 Sq and 2 3 Po states are almost degenerate and, therefore, the PNC ef- 
fect is strongly enhanced. Numerical results are presented for heliumlike gadolinium and thorium, which 
seem most promising candidates for such kind of experiments. In both Gd and Th cases the photon en- 
ergy required will be anticipated with a high-energy laser built at GSI. Alternatively, it can be gained with 
ultraviolet lasers utilizing relativistic Doppler tuning at FAIR facilities in Darmstadt. 

PACS numbers: 1 1.30.Er, 34.80.Lx 
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I. INTRODUCTION 



Measurements of parity nonconservation (PNC) effects with heavy few-electron ions can pro- 
vide new opportunities for tests of the Standard Model at low-energy regime. This is mainly due 
to the fact that, in contrast to neutral atoms (see Refs. 1IH1Q0), in highly charged ions the electron- 
correlation effects, being suppressed by a factor 1/Z (Z is the nuclear charge number), can be 
accounted for by perturbation theory to a very high accuracy. The simple atomic structure of such 
ions allows also one to calculate the QED contributions to the required accuracy. 

PNC experiments with highly charged ions were first discussed in Ref. Ill ill . There it was 
proposed to use close opposite-parity levels 2 1 S and 2 3 Pi in He-like ions for Z ps 6 and Z m 
29, where the PNC effect is strongly enhanced. Later, various scenarios for PNC experiments 
with heavy H- and He-like ions were considered in a number of papers Ill2l-l2lll . In particular, in 
Ref. [13] it was proposed to study the induced 2 3 S' 1 - 2 1 S transition in He-like ions with Z rj 6 
in the presence of electric and magnetic fields. Possibilities to investigate PNC effects in H-like 
ions at high-energy ion storage rings utilizing relativistic Doppler tuning and laser cooling were 



considered in Ref. [17]. Most of the works ||12| . 1141 . 1 161 1 181 . 12 ill exploited, however, the near- 
degeneracy of the 2 1 S and 2 3 P states in He-like ions at Z f» 64 and Z ps 90. For overviews of 
the schemes suggested we refer to Refs. iHT 



210. 



In the present paper, we evaluate the PNC effect on the laser-induced 2 3 S' 1 - 2 1 S transition in 
heavy heliumlike ions nearby Z = 64 (transition energy of about 114 eV) and Z = 90 (transition 
energy of about 240 eV), where the PNC effect is strongly enhanced. Such experiments seem 
to be feasible in near future in view of recent developments in high-energy lasers for heavy ion 



experiments (PHELIX project) li22l |23fl. As an alternative, one may consider employment of 



relativistic Doppler tuning at FAIR facilities in Darmstadt ll24l.l25ll. With ion energies up to 10.7 
GeV/u, as anticipated at the FAIR facilities, the Doppler effect can be utilized for tuning ultraviolet 
laser light with photon energies in the range from 4 to 10 eV in resonance with the transition 
energies under consideration. 

The paper is organized as follows. In sectionHH the basic formulas for the 2 3 Si- 2 1 S , transition 
amplitude are presented. The admixture of the opposite-parity states 2 1 S and 2 3 P is taken into 
account and, as a result, the PNC correction to the cross section is derived. It is shown, that 
accounting for the first-order interelectronic-interaction and QED corrections in the velocity gauge 
can be easily done within the zeroth-order approximation in the length gauge. In section [Till 
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numerical results for the PNC correction in heliumlike gadolinium and thorium are presented and 
possible scenarios for experiments are discussed. 

Relativistic units (h = c = 1) and the Heaviside charge unit [a = e 2 /(An), e < 0] are used 
throughout the paper. 



H. BASIC FORMULAS 



We consider the absorption of a photon with energy to m E 2 i s — E 2 3 Sl and circular polar- 
ization A = ±1 by a heavy heliumlike ion being initially prepared in the 2 3 Si state. If the weak 
electron-nucleus interaction is ignored, the absorption cross section is completely determined by 
the magnetic-dipole transition amplitude. For such a transition the interelectronic-interaction ef- 
fects are suppressed by a factor 1/Z and, to zeroth order, we assume that the electrons interact 
only with the Coulomb field of the nucleus. Then, the wave functions of the initial (2 3 Si) and the 
final (2 1 S'o) state are given by 

«J M (Xl, X 2 ) = -j= ^ C jXuj2 m2 W>31 mi ( X l)^j 2 m 2 (x 2 ) - ^ h mi (x 2 )^ a m 2 (x X )] , (1) 

v mi,m2 

where iftj imi (x) is the one-electron Is wave function, ^ m2 (x) is the one-electron 2s wave func- 
tion, and Cj^ ni j 2m2 is the Clebsch-Gordan coefficient. In what follows, we assume that the laser 
spectral width and the width due to a finite ion-laser interaction time can be neglected. If, for a 
moment, we further neglect the width of the initial state, the cross section in the resonant approxi- 



mation is given by (see, e.g., Refs. [26-28]) 



, 3 r t |(6|[fl(l) + fl(2)]|q)| 2 

a=i2w) (/s, + u -£W» + r>/4- (2) 

Here \a) = |2 3 S'i) and |6) = |2 1 5 , ) are the initial and final states, respectively, T b is the width 
of the final state, and R(i) is the transition operator acting on variables of the ith electron. In the 
transverse gauge, R = —eot ■ A, where 

= eexpj^x) 

v/2w(27r) 3 

is the wave function of the absorbed photon and a. is the vector incorporating the Dirac matrices. 
In order to account for the width of the initial state a in Eq. ([2]), we simply replace E a — )■ E and 

W( ' h / M (g-Sl1/4 ^- W 



In a more rigorous approach, one should consider the preparation of the state a as a part of the 



whole process |27|]. With the substitution ©, we get 

n - frlF T b T a \(b\[R(l) + R(2)]\a)\* 

~ w J dtj 2^\{E + u - E b f + ri/4][(E - E a y + ri/4] ■ ° } 

Integrating over E, we obtain 

- - i 2 ') ! Md^krwi l(4||s(1) + ™ r (6) 

In the resonance case, u = E b — E a , we have 

r a + r& 

Finally, averaging over the angular momentum projection of the initial state, we obtain 

\(b\[R(l) + R(2)]\a}\" 

M a 



(8) 



In what follows, due to smallness of the transition energy, we can write 



R = -ea ■ A = -e(a ■ e) 6X p( ' k ' x) « -e ' e) (1 +ik ■ x) . (9) 

v/2w(27r) 3 V 2w ( 27r ) 3 

For the transition J a = 1 — > J b = we can restrict to the dipole approximation and, therefore, 
represent the transition operator as the sum 

R = R C + R m , (10) 

where 

R c = -e ; J (11) 

is the electric-dipole transition operator in the velocity gauge, 

fl m = i ([6 / Xk] - M) (12) 

is the magnetic-dipole transition operator, and ^ = (e/2)[x x a] is the operator of the magnetic 
moment of electron. If we neglect the weak interaction, the 2 3 S' 1 - 2 1 S transition amplitude is the 
pure magnetic-dipole one. Then, evaluating the matrix elements in Eq. ([8]), we obtain 

vf S ^ = " \(2sM\2s) - (lsM\ls)\\ (13) 

9 J- 23Si + 1 2 l S 
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where (ns\ \ns) is the reduced matrix element of the magnetic-dipole-moment operator and the 
subscript "0" stays for the zeroth-order approximation. 

To account for the weak interaction we have to first modify the wave function of the 2 1 S Q state 

o 



due to the admixture of the 2 3 P n state: 



l^^l^ + MiMMl^. (14) 

Here 

Hw = -(G F /V%)QwPN{r) lb (15) 



is spin-independent part of the effective nuclear weak-interaction Hamiltonian Q29Q . Gp denotes 
the Fermi constant, Qw ~ — N + Z(l — 4sm 2 9\y) is the weak charge of the nucleus (which is 
related to the Weinberg angle 9 W ), 75 is the Dirac matrix, and pn is the effective nuclear weak- 
charge density normalized to unity. A simple evaluation of the weak-interaction matrix element 
yields 

(2 3 P \[H W (1) + H W {2))\2 1 S ) = (2 Pl/2 \H w \2s) 

. Gp 



POO 

Q w drr 2 p N (r)[g 2p f 2s - h Pl/2 92s] , (16) 
Jo 



2^2 

where the large and small radial components of the Dirac wave function are defined by 

/ g nK (r)n Km (n) \ 
VWm(r) = (17) 
\ z/ nK (r)fi_ Km (n) I 

and k = (—iy +l+1 ^ 2 (j + 1/2) is the Dirac quantum number. Then formula (fT4l) can be written as 

|2%)^|2 1 5 ) + ze|2 3 Po), (18) 

where 

5 = ih^-E^ I *^W[ 9 *,„/ 2 .-W 2 J. (19) 

The admixture of the 2 3 P state enables the 2 3 Si — 2 3 P transition, which is determined by the 
electric-dipole amplitude. Since the electric-dipole transition operator depends on the gauge em- 
ployed, the results may differ in the different gauges, if the calculations are restricted to a given 
approximation. The difference can be especially large for a transition between the states hav- 



ing the same (or close) zeroth-order energies, as in the case under consideration. In Ref. OOh 



it was shown that using the length gauge in the calculation of the zeroth-order 2s — 2pi/ 2 tran- 
sition amplitude in H-like ions is equivalent to accounting for the one-loop QED corrections in 
the velocity-gauge calculation, provided the transition energy in the length-gauge calculation in- 
cludes the corresponding corrections. Let us show that accounting for the one-photon exchange 
and one-loop QED corrections to the 2 3 S' 1 — 2 3 P transition amplitude in the velocity gauge can be 
performed equivalently within the zeroth-order approximation in the length gauge by employing 
the transition energy which includes the related corrections. 

With this in mind, we consider first the evaluation of the 2 3 S\ — 2 3 P transition amplitude in 
the velocity gauge to zeroth and first orders in 1/Z. The corresponding diagrams are presented 
in Figs. 1 and 2, respectively. Formal expressions for these diagrams can be derived using the 
two-time Green function method [|27ll . Such a derivation was considered in detail in Ref. [31]. To 
simplify the analysis, we consider the matrix elements of the electric-dipole transition operator 
between the one-determinant wave functions, 

Uo(xi,x 2 ) = -^=^(-l) p ^ Pai (x 1 )V'Pa 2 (x2), (20) 

«&(xi,x 2 ) = -^=^(-l) p ^p fel (xi)^p 62 (x 2 ) , (21) 

where it is assumed that ai = b± — Is, a 2 — 2s, b 2 = 2py 2 , P is the permutation operator, and 
( — l) p is the sign of the permutation. Then, to zeroth order we obtain for the transition amplitude 

r (°) = -<6|[i2e(l)+i2e(2)]|a> = -<&i|P e (l)|ai)(5 a2 62 - (b 2 \R c (2)\a 2 )5 aibl 

= -(2 Pl/2 \R c \2s) . (22) 

Employing the identity 

a = i[H,r], (23) 
where H is the one-electron Dirac hamiltonian, one obtains 

r<°) = i 7 =4= (2 ft/2 |(e • r)|2 S )( £2 - e 2s ) , (24) 

where e 2s and £2 Pl/2 are the one-electron Dirac energies of the 2s and 2p%/ 2 states, respectively. 
In particular, it follows that for the pure Coulomb field (e 2s = s 2pi/2 ) in the velocity gauge the 
zeroth-order 2 3 Si — 2 3 P transition amplitude is equal to zero. 
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The interelectronic-interaction corrections j_defined by the diagrams depicted in Fig. 2, consist 



of irreducible and reducible contributions [27, 



3 ill . Since, to a good accuracy, these corrections can 



be treated with the pure Coulomb field of the nucleus, in what follows, we restrict our consideration 
to this approximation. Then, according to Ref. Bill we find that for the 2 3 5*i — 2 3 P transition the 
reducible contribution vanishes. As for the irreducible contribution, it can be expressed as the sum 



L31] 



where 



.(a) 



_ (a) , (6) 
7~irr T; rr T Tj rr , 



(25) 



1 



£(-l) P { £ (Ph\(e-cx)\n) 



x- 



(nPb 2 \I{e Pb2 - Ea^la^) 



(») 



+ ( pfo 2|(e-a)|n)^ (PMI-^p&i - e ai )|aia 2 )} , (26) 



irr 



2w(2 



7T 



£(-!) p { E (w^i/^ - 

P n 



£a 2 )l na 2) 



X • 



(n\(e ■ a)\ai) 



+ (PbiPb 2 \I(ep bl - e ai )\am) 



E (0) -e -e 



(n|(e-a)|a 2 )}. (27) 



Here I(u) = e 2 a^a u D 



P> pCT (w,x- y) = -g pc 



dk exp (zk • (x — y)) 
(2tt) 3 w 2 - k 2 + iO 



(28) 



is the photon propagator in the Feynman gauge, a p = 7°7 P = (1, a), Pi -* = e ai + e a2 , and 



E^ = e 6l + £ b2 . Taking into account that = E^ and using the identity (1231) . we get 

(^Pfe 2 — £ a 2 )l a l a 2/ 

■e 01 )|oia a )}, (29) 



(a) 



171 



E(-!) P { E (^i|(e-r)|n)(nP6 2 |/( 



V / 2cJ(2 

+ ^ (P6 2 |(e-r)|n)(P6 1 n|J(£ P6l 
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-(b) e 



= /^==^(-l) p { ]T (P6 1 P6 2 |J(£p 62 -£ a2 )|na 2 )(n|(e-r)|a 1 ) 



p 



/ 2cu(2 7 

+ ^ (PfeiPfeal^ep^-eaJlamXnKeT)!^)}. (30) 

n 

With the aid of the completeness condition 

J2\n)(n\ = l, (31) 



we find for the sum of the expressions (1291 ) and (|30 



r irr = ^ . % = (2p 1/2 |(e ■ r)|2s)(AP b - AP a ) , (32) 
a/2w(27 



where 



AP a = 53(-l) p (Pa 1 Pa 2 |/(e Pai - e Bl )|aia 2 ) , (33) 
p 

AE b = 53(-l) p (P6 1 P6 2 |/(£ P6l - ejJI&x&a) (34) 
p 

are the first-order interelectronic-interaction corrections to the initial and final states, respectively. 
The same relation holds if one includes the one-loop QED corrections. The corresponding proof, 
which was given first in Ref. [3(3], is presented in the Appendix. Summing up the zeroth- and 
first-order contributions yields 



T = <vm7 {2piM(e ' t)]23HEi " E - ] = l 7m? {2piM[e ' d)|2s> • <35) 



where d = er is the operator of electric-dipole moment, E a and E b are the total binding energies of 
the initial and final states, respectively. It is evident that similar equations can be derived involving 
two-electron wave functions (OQ). Consequently, in what follows, we will use the electric-dipole 
transition operator in the length gauge 

^ = . l J£t£L. (36) 



Substituting the two-electron wave function (TT8T) into Eq. © and performing the calculation, we 
obtain for the PNC contribution to the cross section 

4nc" 2150) = " 2\£(2p 1/2 \\d\\2s)({2s\\»\\2s) - (ls|M|ls>) , (37) 



where (2pi/ 2 ||d||2s) is the reduced matrix element of the electric-dipole-moment operator and 
A = ±1 is the photon polarization. Integrating over the angular variables in the reduced matrix 
elements yields 

poo 

(ns\\n\\ns) = -2ey^2/3 / dr r 3 g ns (r) f ns (r) , (38) 

Jo 

POO 

(npi/ 2 ||d||ras) = -e^/2/3 / dr r 3 (g npi/2 (r)g ns (r) + f npi/2 (r)f ns (r)) . (39) 

Jo 

These integrals are easily evaluated employing the virial relations for the Dirac equation (see, e.g., 



Refs. 132 



34(|). For the case of interest here, one derives 



drr 3 (g 2s (r)f 2s (r)-g ls (r)f ls (r)) = 7- y/(l + i)/2, 



drr 3 (g 2p (r)g 2s (r) + f 2piJ2 {r)f 2s {r)) 



3(1 + 7)^^+2^ 
2aZ 



where 7 = ■sjl — (aZ) 2 . Substituting these expressions into Eqs. (fT3l) . (1371) leads to 



a 



o. 



(2 3 5i^2 1 S ) 



(2 3 5i^2 1 5 ) 



nauo 



27T 2 3. Sl + r 2 i5 



(V2(l + 7)-2 7 ) : 



cr, 



(2 3 Sn-2 1 S ) 1 (2 3 5i->2 1 5 ) 







+ <J 



PNC 



(2 3 5i^2 1 S ) 



where 



2e 



(2P1/2 1 


\d\\2s) 


(2s\ 




\2s) - 


-( 






|ls> 



-2e- 



3(l + 7 )v^ + ^ 



(40) 
(41) 

(42) 
(43) 

(44) 



'aZ(v/2(l + 7)-2 7 ) 

is a parameter which characterizes the relative value of the PNC effect. The second term in the 
right-hand side of Eq. (1431) represents the PNC contribution, which changes the sign under the 
replacement A — ?- — A. The PNC parameter can also be represented as 



-26 



T2 3 5i + T2 3 P O"o 



(2 3 SW2 3 P ) 



\ r 23Sl + rv So a fs^ So) 



(45) 



where a, 



(2 3 Si^2 3 P ) 



is the cross section of the resonant absorption into the 2 3 P state and r 2 3p is 



the total width of this state. 



in. RESULTS AND DISCUSSION 



The formulas (l4~2l)-(l44l allow one to evaluate the cross section and the corresponding PNC 
effect. The most promising situation for observing the PNC effect occurs in cases where the levels 
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2 1 Sq and 2 3 P are almost degenerate. According to the most to-date elaborated calculations 05 - 



3811 (see also the related table in Ref. [21]) such cases are gadolinium (Z = 64) and thorium (Z 



90), where the levels 2 1 S and 2 3 P are near to cross. In case of Gd the energy interval amounts to 
-0.023(74) eV, while in case of Th it is 0.44(40) eV B. Since the uncertainties are comparable 
to the energy differences, to estimate the PNC effect we take the values 0.074 eV and 0.44 eV 
for the 2 3 P - 2 1 So energy difference in Gd and Th, respectively. The widths of the 2 3 Si and 
2 1 So states, which enter formula (|42l) . are mainly defined by the one-photon Ml and two-photon 
E1E1 transitions, respectively. We evaluate the decay rate of the Ml -transition 2 3 Si — > l 1 ^ 
employing the transition energy taken from Ref. 0511 . The interelectronic-interaction corrections 
to the transition amplitude are calculated to first order in 1/Z within a systematic QED approach 
(see for details Ref. [31]). As the result, we obtain decay rates w^ 31 ^ 1 = 2.301 x 10 12 s" 1 for 
Gd and w^u^ 1 ^ = 9.470 x 10 13 s _1 for Th. These values are in a fair agreement with those 
from Refs. HflB- The two-photon decays 2 3 S 1 ->• 1% and 2% 1% are calculated in 
the length gauge with the transition energies taken from Ref. |p5|l . The interelectronic-interaction 
effects are approximately accounted for by means of a Kohn-Sham potential. The calculated 
transition rates are w ( ^ Sl ^ 11So) = 8.74 x 10 8 s" 1 , w ( £ So ^ 11So) = 9.04 x 10 11 s" 1 in case of 
Gd and w^ Sl ^ 11So) = 2.07 x 10 10 s~\ w { £ S °^ 11So) = 6.25 x 10 12 s" 1 in case of Th. These 
values together with the E1E1 channel include also higher multipole contributions, such as M1M1 
etc. In case of Th the dominant E1E1 decay channel yields ^eiei"^ ^ = 1-62 x 10 10 s _1 and 
u, EiEi°~ i ' 1 = 6.25 x 10 12 s -1 . It is worth noticing that for the 2 3 5i state the higher multipoles 
contribute up to 20% to the total two-photon decay rate. Comparing the E1E1 decay rates with 
the results of Ref. |4l|], we find an excellent agreement for the 2 3 Si state and a slight deviation 
for the 2 1 S'o state, which is mainly due to employing the more accurate transition energies in our 
calculations. Finally, the total widths are T 2 ss 1 = 1.515 meV, r 2 is = 0.595 meV in case of Gd, 
and r 2 3£ 1 = 62.35 meV, r 2 i5 = 4.11 meV in case of Th. The results of the calculations of the 
PNC effect by formulas (O-dS]) for Gd and Th are presented in Table I. 

As one can see from the table, in both Gd and Th cases the PNC effect amounts to about 0.05%, 
which is a rather large value for parity-violation experiments. What is more, one may expect that 
the PNC effect can be further increased, at least, by an order of magnitude by choosing proper 
isotopes, provided the 2 1 S — 2 3 P n energy difference is known to a higher accuracy. With the 
current experimental techniques [42], accurate measurements of the difference considered seem 
feasible. 
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Because of a large transition energy (> 100 eV), until recently the experimental scenarios with 
the laser-induced 2 3 Si — 2 1 So transition in heavy He-like ions were far from being possible. How- 
ever, the situation has changed in view of the very significant progress in X-ray laser development. 
Such lasers will be available in the near future with a high repetition rate [43]. Already now, 
there is a first X-ray laser available at the heavy ion facility GSI (PHELIX facility) where photon 



230. As an alternative scenario, the excitation 



energies of up to 200 eV have been reached [22, 
energy can be obtained by counter-propagating the ultraviolet laser beam with the photon energy 
in the range from 4 to 10 eV and the He-like ion beam with the energy up to 10.7 GeV/u, which 
will be available at the FAIR facility in Darmstadt 11241 12511 . The population of the 2 1 S level can 
be measured by observing the 2E1 decay to the ground state. In the second scenario, due to the 
strong Lorentz boosting, the decay photons are emitted at the forward direction, that considerably 
simplifies their detection. 

The next problem to be addressed is the preparation of ions in the 2 3 S\ state that is required in 



both scenarios considered. As follows from the study presented in Ref. 11441 14511 . in collisions with 
gas atoms one can produce selectively both the 2 1 S state and the 2 3 S± one |46|]. However, it would 
be of great importance to populate exclusively only the 2 3 Si state. The only way to accomplish 
this is to form first the doubly excited (2s2pi/ 2 )o state via dielectronic recombination of an electron 
with a H-like ion. Since the main decay channel of the (2s2pi/ 2 ) state is the transition into the 
2 3 Si state, this enables selective production of ions in the 2 3 S' 1 state. 

The PNC effect is to be measured by counting the intensity difference in the 2E1 decay of the 
2 1 So state for polarizations A = ±1. The background emission can be separated by switching 
off the laser light. Changing the photon energy allows one to eliminate the interference with a 
non-resonant transition via the 2 3 P state, which could also be evaluated to a good accuracy if 
necessary. Moreover, since the 2E1 emission can be measured relative to the intensity of the Ml 
X-ray line (decay of the 2 3 S' 1 state), such an experimental scenario appears to be quite realistic. 



IV. CONCLUSION 

In this paper we have studied the PNC effect with laser-induced 2 3 S' 1 - 2 1 S transition in heavy 
heliumlike ions. A simple analytical formula for the photon-absorption cross section derived en- 
ables easy evaluation of the PNC effect for ions nearby Z = 64 and Z = 90, where the effect 
is strongly enhanced due to near-crossing of the opposite-parity 2 1 Sq and 2 3 P levels. The cal- 
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culations performed showed that the effect can amount to about 0.05% and even bigger for the 
ions of interest. Prospects for the corresponding PNC experiments have been discussed. It is 
found that the desired photon energy can be achieved either by X-ray lasers that are presently 
getting developed at GSI (PHELIX project) as well as at the Helmholtz-Institute in Jena Il47ll or 
by counter-propagating the ultraviolet laser beam and the He-like ion beam at the FAIR facility in 
Darmstadt . 
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Appendix: QED corrections to the transition amplitude 



The one-electron QED corrections to the 2 3 Si — 2 3 P transition amplitude are determined by 
the corresponding contributions to the 2s — 2pi/ 2 amplitude in one-electron ions as defined by 
the diagrams shown in Fig. 3. Formal expressions for these corrections are almost the same as 



for the corresponding corrections to the emission amplitude 112711 . Let us consider the one-loop 
self-energy correction. According to formulas provided in Ref. 112711 . it is given by the sum of the 
irreducible, reducible, and vertex contributions. For the electron interacting with pure Coulomb 
field together with the dipole approximation exp (ik ■ x) — > 1, the reducible contribution vanishes. 
The irreducible contribution is given by 

rif } = -{2p 1/2 \R\£ 2s ) - {& Pm \R\2s) 

= [<2pi/2|(e • oOlGu) + (£apj(e ■ <x)\2s)) , (A.l) 



y / 2u(27l) 3 ' 



where 



a \n)(n\E(s a )\a) 



and 



fc,» = E _ — ■ < A - 2 > 



(A.3) 



a\ne)\b) = if ' duY H^^'f . (A.4) 



2tt7-oo V e-w-e n (l-iO) 
By means of the identity (|23l) and the completeness relation (I3TI) . we obtain 

7^f r E) = i ^__ {(2p 1/2 |(e • r)\2s)({2 Pl/2 \Z(e 2pi/2 )\2 Pl/2 ) - (2 S |S(e 2s )|2 S )) 

+(2p 1/2 |[(e-r),E(£ 2s )]|2 S )}. (A.5) 

For the vertex contribution one derives 



ni,n 2 



1 e(e • a) 1 



e 2pi/2 - w - £ ni (l - iO) v /2w(27r) 3 £ 2s - a; - £ n2 (1 - iO) 
x(n 2 \apexp (— ik • x)|2s) . (A. 6) 
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Transforming 



£2 Pl/2 -u}-e ni (l- iO) e 2s -u- £ n2 (l - *0) 



£ ni -^n 2 ^2 Pl/2 -w-e ni (l -z'O) £ 2s -w-e na (l -z'O), 



(A.7) 



we get 

r (SE) 



(ni|(e • a)|n 2 ) = (e • r)]|n 2 ) = z'(£ m - £n 2 )("i|(e ■ r)\n 2 ) 



(A.8) 



|3 2tt 



da; 



dk 



xz 



(2tt) 3 a; 2 - k 2 + zO 
1 



(2pi/ 2 |a p exp (ik-y)|ni) 



ni,n 2 



£ 2 Pl/2 - w - £ m (1 - z'O) e 2s -uj- e n2 (l - z'O) 
x (n 2 |a p exp (— zk ■ x)|2s) 



(ni|(e • r)|n 2 ) 



v/2a;(27r 



e „ i i°° , /" dk 
' da; 



|3 2tt 

1 



(2vr) 3 u- 



_ k 2 + i0 {l^( 2 ^ 1 /2l QiPex P (* k • y)lm> 



(rii | (e • x)a p exp (— z'k • x)|2s) 



£ 2 Pl/2 -w-e m (l-iO) 
- (^i|( e • r )l^2)(^2|a P exp (-ik • x)|2s) 

"2 



I (2p 1/2 1 (e ■ y)a p exp (ik • y) \n 2 ) 

(2pi /2 |a p exp (ik • y)|m)(ni|(e • r)|n 2 ) 



x 



e 2s - uj - e n2 (l - z'O) 



— (n 2 |a p exp (-ik ■ x) |2s) | 



y/2u(2n) 



:<2p 1/2 |[(6-r),E(e 2a )]|2s). 



(A.9) 



The sum of both irreducible and vertex contributions yields [13011 



(SE) 



tot 



v/2o;(2 



7T 



(2p 1/2 |(e • r)|2 S )((2p 1/2 |S( £2pi/2 )|2 Pl/2 ) - (2 S |E( £2s )|2s)) . (A.10) 



A similar equation can be derived for the vacuum-polarization contribution. 
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